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Abstract. We consider a modification of the fluid flow model for a swirling vortex developed by J. Serrin 
[23| , where velocity decreases as the reciprocal of the distance from the vortex axis. Recent studies, based 
on radar data of selected severe weather events [2| |28||29| , indicate that the angular momentum in a tornado 
may not be constant with the radius, and thus suggest a different scaling of the velocity /radial distance 
dependence. 

Motivated by this suggestion, we consider Serrin's approach with the assumption that the velocity de- 
creases as the reciprocal of the distance from the vortex axis to the power b with a general 6 > 0. This leads 
to a boundary- value problem for a system of nonlinear differential equations. We analyze this problem for 
particular cases, both with nonzero and zero viscosity, discuss the question of existence of solutions, and use 
numerical techniques to describe those solutions that we cannot obtain analytically. 



1. Introduction 

In 1972, J. Serrin 23 , following the works of Long [l8 19 and Gol'dshtik [l2], discovered a special class 
of swirling vortex solutions to the Navier-Stokes equations, where velocity decreases as the reciprocal of the 
radial distance, r, from the vortex axis. Depending on kinematic viscosity and the value of a "pressure" 
parameter, he described the cases of downdraft core with radial outflow, downdraft core with a compen- 
sating radial inflow, and updraft core with radial inflow (single-cell vortex). Search for axisymmetric flow 
solutions has continued through the last few decades because of their importance in modeling a wide range 
of phenomena. Relevant reviews are presented in [25(|26) . 

However, recent radar studies indicate (see, e.g., [2]) that tornado-related vorticity flelds might have 
fractal nature with respect to the grid size; more specifically, natural log of the vorticity, C, and natural 
log of the grid spacing seem to have a linear relationship, with a constant negative ratio. Larger absolute 
values of this ratio correspond to stronger storms. In some cases of tornadic storms, the ratio is found close 
to —1.6. For tornadic mesocyclones, this suggests a power law of the form ^ (x . Cai |2] observed that 
the exponent can be thought of as measuring a fractal dimension associated with the vortex. In a related 
work |6 , we conjecture that stretching of the vortex lines could lead to a fractalization of the vortex and, 
possibly, to a vortex breakdown. Additional discussions of vortex breakdown can be found in p]|8pp5}|T7pO] . 
We comment on how our results relate to vortex breakdown in the conclusions section. The possibility of 
fractalization in the process of vortex stretching was pointed out by Chorin [4|. 

For high-resolution mobile Doppler radar data, there has been some attempt to interpret Cai's power law 
as a power law for the drop-off of the velocity as a function of the radial distance from the vortex axis |28^|29j . 



In these papers, devoted to analyzing data associated to strong or violent tornados, an atte mpt is made to 
calculate the value of the exponent in the "velocity power law" v oc r^. Wurman and Alexander [29' calculated 
the exponent b from radar data obtained in an intercept of the May 31, 1998 Spencer, South Dakota, tornado 
and obtained the value b = —0.67 for the velocity field away from the core-flow region, in which data indicated 
a solid-body rotation. The tornado was rated EF4. These values were calculated from the data taken at one 
instant during the tornado's existence. In the case of the June 2, 1995 Dimmit, Texas, tornado, Wurman and 
Gill [28] observed exponent values in the range —0.5 to —0.7, concluding that "it implies that the angular 
momentum in the tornado was not constant with radius, but decays toward the center." 
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We therefore find it natural to ask whether there are Serrin-type similarity solutions to the Navier-Stokes 
equations with the velocity field proportional to l/r**, where & > and b 1, with the most interesting 
case being < 6 < 1. This work attempts to answer this question. We show that under the assumption 
of constant nonzero viscosity and suitable assumptions on the form of the velocity field, similar to that 
in [23j, the Navier-Stokes equations do not admit any nontrivial solutions except when 6=1, the case 
studied by Serrin. However, in the relaxed case of zero viscosity, the Euler equations always admit a simple, 
purely rotational solution with azimuthal velocity of the form C/r^. In addition, when 6=1, another set 
of nontrivial solutions is found analytically. When 6 > 2, we show that no other solutions exist. The most 
intriguing cases are when < 6 < 1 and 1 < 6 < 2, for which we have not found explicit solutions; for the 
former case we present numerically computed solutions supporting our claim of existence of solutions, while 
for the latter case we show that any solution would have to be unstable in the sense of Rayleigh's circulation 
criterion 

The paper is organized as follows. In section [2] we describe the basic geometry of the problem, the 
governing equations, and the form of solutions we are interested in finding. In section [3] we analyze the 
Navier-Stokes equations in the case of constant nonzero viscosity. In section |4] we focus on the case of zero 
viscosity, governed by the Euler equations. Finally, in section [5] we discuss the conclusions and implications 
of our findings for tornadogenesis. The appendix contains some auxiliary equations needed for our work that 
would unnecessarily clutter the presentation in the paper. 

2. Governing equations, basic geometry, and modified Serrin's variables 

In this section, we discuss the relevant equations, introduce a change of variables in the spirit of |23|, and 
also introduce a special form of solutions we seek, which eventually allows us to reformulate the problem in 
terms of ordinary differential equations. Finally, we discuss the continuity equation and its implications in 
terms of boundary conditions. 

2.1. Governing equations. The equations governing fluid flow are the Navier-Stokes equations, 

P^ = P[^ + i^-^)^) = + (-^ + P)^i^ • v) + Av, (2.1) 



Dt \dt 

where v, P, and p are velocity, pressure, and density fields, respectively, and fi and A are dynamic viscosity 
coefficients. 

The conservation of mass, or continuity, equation is 

^+pV-v= |^+vVp + pV-v = 0. (2.2) 

We will consider the case of an incompressible flow, so that Dp/Dt = V • v = 0. We will seek steady-state 
solutions, i.e., those that satisfy dv/dt — dp/dt — 0, which in turn yields Vp = 0, and therefore p is constant 
and the flow is homogeneous. As a consequence, steady-state solutions for an incompressible flow satisfy the 



following simplified versions of (2.1) and (2.2) 



(vV)v = -Vp-^i^Av (2.3) 

and 

V • V = 0, (2.4) 
where p = P/p is a, (scaled) pressure field and v = p/ pis a, (constant) kinematic viscosity. 

2.2. Spherical coordinate system and components of the velocity field. Following Serrin f23^, we 
use the (right) spherical coordinates {R,a,6), where R is radial distance from the origin, a is the angle 
between the radius vector and the positive z-axis, and 9 is the meridian angle about the z-axis. The positive 
z-axis then corresponds to a = and the boundary (ground) plane to a = 7r/2. We are interested in solutions 
in the upper half space, z > 0, which corresponds to i? > and < < 7''/2 in our coordinate system. 

We denote the components of the velocity vector v{R,a,6) in the spherical coordinate system by vj^, Va, 
and vg, and write 

v(i?, a, 9) = {vr{R, a, 9), Vo.{R, a, 9), ve{R, a, 9)) . (2.5) 
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We will refer to the individual components as radial (wfl), meridian and azimuthal {vg). The scaled 

pressure field will be denoted by p{R,a,9). The three components of the Navier-Stokes equations (2.3) in 
this coordinate system and in the spherical velocity components are given in the appendix in (6.1|-(6.3l, 
and the continuity equation (2.4 1 is given in (6.4 1. 

We will follow Serrin's approach |23j and consider velocities of the form 



VR{R,a,6) 



G{x) 



F(x) 



where r 



r" r" 
Rsma is the horizontal distance to the z-axis, x 



ve{R, a, 6) 



n{x) 



(2.6) 



cos a, and & > 0. We remark that the case 



studied by Serrin 23 corresponds to 6 = 1. 

Since sin a = v 1 — x"^, we can use the change of variables 

fix) = Fix)il - g{x) = G(a;)(l - x^Y'-'^', 

and rewrite ( |2.6| as 

9i^) lu /l^ _ fi^) 



vif\R, a, ( 



R^ si 



sma 



Va{R, a, 0) 



uj{x) 



R^ sn 



ve{R,a,9) = 



Rb sn 



(2.7) 



When 6=1, the upper-case functions, F, G, and D,, agree with the lower-case functions, /, g, and uj. 



The continuity equation (2.4 1 written in terms of the newly introduced functions is given below in (2.8 I, 



while the Navier-Stokes equations (2.3 1 are discussed in section 2.4 



We also remark that the function 

*(i?,a) = R^^''f{x) = i?2-''/(cosa) 

satisfies 

V* • V = 0, 

and therefore the surfaces ^I* = constant contain the streamlines of the fluid motion. Notice that this is a 
direct generalization of Serrin's ^E* — RF{x), since when & = 1 we have F{x) = f{x). 



2.3. The continuity equation. We now consider the continuity equation (2.4) and its consequences for 
solutions to (2.3). We first observe that direct substitution of (2.6) into the continuity equation (2.4) yields 

(2 - h)G{x) ^ ^Jl-x'^ F'{x) - (1 - b) , ^ Fjx), 



(see (6.8) in the appendix) 



Vr 



{2~b)gix)= ^l-x^f{x). 



(2., 



The prime symbol will denote differentiation with respect to x throughout the paper. From (2.8) we see that 
if 6 7^ 2, then G can be expressed in terms of F and g in terms of /. 

We next derive an integral version of the continuity equation that will lead to boundary conditions 
needed later. Let i?o > and £' C M'^ be the upper half ball bounded below by the horizontal disk 
D = {{R,a,0) : < R < Rq, a ^ tt/2, < 6 < 2t:} and above by the hemisphere S ^ {{R,a,e) : R = 
Ro, 0<a< 7r/2, 0<9 < 2tt), i.e., 

E^{{R,a,e): 0<R<Ro, 0<a< 7r/2, < 6* < 2t:}. 



Applying (2.4) and the divergence theorem, we obtain 

0^ Iff V vdT/= // v-ndA 



dE 



Vn dA 



However, since there can be no source or sink at the ground (a 
F(0) /R'' = for all i? > 0, or 

F(0) = /(O) = 0, 



and thus the integral over the disk D vanishes. Substituting (2.6) into the integral over S, we obtain 

r-/^Gix 



Vr dA 



s 



- Rl sin a da de = 27ri?^"^ 



Vr dA. 

7r/2), we have that Va{R,7T/2,9) = 

(2.9) 

le integral over S, we obtain 

tt/2 



0(2;) (sin a) '""da = 
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or, in terms of x, 



G{x) 



{l-x^f/: 



dx 



dx = 0. 



It follows from (2.8 1 and (2.9) that if 6 7^ 2, then (2.101 can be rewritten as 

lim F{x){l - a;2)(i-'')/2 lim f{x) = 0. 



(2.10) 



(2.11) 



Remark 2.1 (Consequences of the continuity equation). We can summarize the consequences of the conti- 
nuity equation (2.4| as follows. 



(1) Ifb^2, then, using ( |2.10[ ), and , we have 

G{x) 



F{0) = 0, 
and, in particular, F ~ f = Q. 



dx = 0, 



F{x){l 



-1/2 



const, 



(2) Ifb^2, then, using ( |2.11[ ), and (|2^, we have 

F(0) = hm F{x){l - a:2)(i-'')/2 =0, (2 - b)G{x) = y^l-x^ F' (x) - (1 - &) 



Fix), 



or, in terms of the lower-case functions, 
/(O) = lim fix) = 0, 



(2 - 6).g(x) = - a;2 fix). 



Remark 2.2. Note that the special case b = 1 in the previous remark gives rise to i^(0) = 0, Ymix^i Fix) — 0, 
and Gix) = \/l — a;^ F'ix), which is consistent with 



components have the form w(i?, a, ( 



2.4. Simplification of the Navier— Stokes equations. In this section, we consider the Navier-Stokes 
equations (2.3 1 in terms of the velocity expressions (2.6 1 and (2.7 1. We first observe that all of the velocity 

Kia)/R^, so their partial derivatives are of the form 

dv Kia) dv 



dv Kia) 



= 0. 



da ' 89 

The dot symbol will denote differentiation with respect to a throughout the paper. Note that all terms in 
the left-hand sides of (6.1|-(6.3l, arising from the convective term in (2.3l, are of the form C(a)/i?^+^'', 
while all terms in the right-hand sides, arising from the diffusive term in (2.3 1, (i.e., all the terms multiplied 
by the viscosity coefficient, v) are of the form Dia)/R^+^. Therefore, the governing equations (|6.1|-(|6.3| 
have the general form 

1 dp , ^2.13) 



i?l+26 OR i?2+b ' 

C2(a) _ _}^di^ D2ia) 

Rda 
1 dp 

R sin a 89 



Dsia) 



(2.14) 



Csja) ^ 

i?i+2b Rsinad9^'^ R'^+'> 
where the expressions for Ci(a) and Diia) are given, in their various forms, in the appendix. 

Like in |23], we argue that (2.141 yields dp/ 89 independent of 9, so p must be linear in 9. Together with 
periodicity in 9, this implies that 8p/89 = 0, so piR,a,9) =p(_R, a). Consequently, (2.14) reduces to 



C^ia) Dsia) 



(2.15) 



The scaled pressure function has to satisfy (2.12), so by integrating it with respect to R we obtain 

Diia) 



piR,a) 



Cija) 
26i?26 



+ T(a). 



il + b)R^+b 

Substituting this expression into (2.13), we conclude that Tia) — 0, and thus T(q;) s T is a constant and 

Giia) Diia) 



piR,a) = 



26i?2f' 



(1 + b)R^+'' 



T. 



(2.16) 
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In addition, from (2.13) and (2.161 we obtain a compatibility condition for the existence of the pressure 

Ci{a) + 2bC2{a) _ Di{a) + {1 + b)D2{a) 



2&i?i+2f' 

We now have the following equivalence lemma. 



{l + b)K 



'2+6 



(2.17) 



Lemma 2.1. The system of Navier-Stokes equations (2.12 |-(2.14l is equivalent to the system (2.15|-(2.17). 



Proof. First, note that (2.15l-(2.17l follow directly from the Navier-Stokes equations (2.12 1-(2.14). Vice 
versa, iiCi{a) and Di{a) are such that (2.15) and (2.17) are satisfied, then (2.16) gives an expression for the 
scaled pressure so that (2.12) is immediately satisfied, (2.13) is satisfied due to the compatibility equation 
(|2.17|, and \2.1A\ follows immediately from (|2.15|. □ 



Remark 2.3. In what follows, we will focus on finding solutions to (2.15) and (2.17). We will rewrite and 
use them in a slightly different form, 

C3(a) = vB^-^D^ia), (2.18) 

and 

2b 



C\{a) + 2bC2{a} = lyR 



h-l 



Di{a) + {l + b)D2ia) 



(2.19) 



1 + 6 

3. Analysis of the viscous case: u > 

In this section we discuss the existence of solutions to our problem with constant nonzero viscosity. We 
show that nontrivial solutions of the form (2.6) exist only for the case 6=1 discussed by Serrin |23 . We 
start by discussing the boundary conditions and then analyze the various cases that arise for various values 
of 6. 

3.1. Boundary conditions. In the case of nonzero viscosity, the no-slip requirement at the ground implies 

VR{R,7:/2,e) = ve{R,7r/2,0) = for all i? > and < 6* < 2tt, so that, using G(0) = f2(0) = 0. 

The no-sink/source requirement gives Va{R,'n'/2,d) — 0, or (2.9), F{0) — 0. In addition, as a consequence 
of incompressibility, in particular (2.8), we have F'{0) — 0. (Note that as discussed in Remark |2.1[ if 6 = 2, 
then F = and all boundary conditions concerning F are automatically satisfied.) 

Near the vortex axis, we have (2.11) if 6 7^ 2, while if 6 = 2, then F = 0, and there are no physical 
restrictions on the behavior of G and as x — > 1. However, we will assume, similarly as in [23j , that near 
the vortex axis the azimuthal velocity, vg, behaves like C/r^, i.e., we will assume that 

lim n{x) = lim uj{x){l - x^)-'-^-''^/^ ^C^^Q. (3.1) 

x—^l a:— >-l 

Remark 3.1 (Boundary conditions in the viscous case). In the case of constant nonzero viscosity, sought 
solutions F, G, and are subject to the following requirements: 
• Ifb = 2, then 

F = / = 0, G(0) = n{0) = g{0) = w(0) = 0, lim n{x) = G^. 



Ifb^2, then 
F(0) = F'{0) = G(0) = n{0) = 0, 



lim F{x){l 

X—^l 



p2)(i-b)/2 lim ri(:E) = C^. 



/(O) = /'(O) = .9(0) = c^(O) = 0, lim/(x) = 0, lim c^(a;)(l - x^)-^-'')/^ = G^. 

a;— >1 a;— >1 

Remark 3.2. Note that the boundary conditions studied by Serrin are consistent with ours when 6=1, 
since then lima-^i F{x){l - a;2)(i-'')/2 = lim^^i F{x) = 0. 



3.2. Governing equations. The governing equations are (2.18) and (2.19), together with the continuity 
equations (2.8). We will need to distinguish between the special case 6 = 1 studied by Serrin and the 
remaining cases when 6^1. 
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3.2.1. Case 6=1. In this case, (2.181 and (2.191 imply 



Ca{a) = vDz(a), Ci{a) + 2C2{a) ^ u{bi{a) + 2D2{a)). 



Using (6.9l-(6.12), these equations yield 



1/(1 - - ^vxF"'{x) + F{x)F"'{x) + ZF'{x)F"{x) = 

1 — x^ 

v{l - x'^)Q."{x) + F{x)Vl'{x) = 0. 



(3.2) 



We note that system (3.2 1 is identical to system (5) in ^23j and is analyzed there. In what follows, we focus 
on the case with 6 7^ 1. 



3.2.2. Case b^l. In this case, the only way to satisfy (2.181 and (2.19) for all i? > is to satisfy 

C3(a) = 0, £'3(0) =0, CiCa) + 2feC2 = 0, I)i(a) + (1 + 6)Z?2 = 0. (3.3) 
The relevant expressions for the quantities in (|3.3[) are given in the appendix and we will recall them as 



(2 - h)g{x) = Vl - x^f'ix). (3.4) 



needed. The last equation that needs to be satisfied is (2.8), the consequence of the continuity equation, 
restated here in both forms for completeness, 

{2-h)G{:x) = ^/l-x'^F'{x)-{l-b) F{x) or 

V 1 — 



3.3. Case 6 = 2 (no solutions). We first address the special case with 6 = 2. In this case, F = f = Q and 
(3.4) provides no information. Using (6.5l-(6.7l, the first three equations in (3.3) reduce to 

(3.5) 
(3.6) 

(3.7) 



Gix)n{x) = 0, 

(1 - x^ fn"{x) + 2x{i - x^)n'{x) + 3n{x) = o, 

8a; 



1 



■ G\x) + 2{G\x)y + (n^x))' = 0. 



Because of the boundary condition (3.1 1, 7^ on some interval {xq, 1) by continuity. Equation ( |3.5[ ) then 
implies that G = in (a;o,l), and (3.7) reduces to (ri^)' = in (a;o,l). Thus fi = C^i in (a;o,l). This, in 
turn, implies that Q = C^ and G = in (0, 1). However, neither ( |3.6[ ), nor the boundary condition 17(0) = 
are then satisfied, so no solution exists when 6 = 2. 



3.4. Case 6 7^ 1,2 (no solutions). In this case, (3.4) can be substituted into (3.3) to yield the set of 
equations given in (6.13 1-(6.16 1 in the appendix. 

In order to conclude that no solutions exist, it suffices to analyze ( 6.14[ ), which reads 

(1 - xYn"{x) - 2(1 - 6)a;(l - x^)n'{x) - (1 - b^)n{x) = 0. 

When equipped with the initial conditions 51(0) = and fl'{0) = G, its solution is 

1-6 63 
2^' 2' 2' 

where 2-7^i is the Gaussian hypergeometric function given by (see [21] ) 
and {x)n with n €E N U {0} is the Pochhammer symbol given by 

'l if 71 = 



n{x)^Cx{l^x^)^'-'y\F,{^-^\-l;x^^ 



(3.8) 



{x)n = 



x{x + I) . . .{x + n - I) if rt > 0. 



Since the case G = in (3.8) would yield £7 = 0, which does not satisfy the boundary condition (3.1), we 
only consider the case with G 7^ and show that (3.1 1 cannot be satisfied for any 6 > 0. 
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We first have 21 that 



Hm 



1-663 . 



r(i)r(i) 
r(¥)rm 



2r(^)r(^) 



Since 1/T(z) is an entire function vanishing only for z = 0, —1, —2, . . . , the limit above is finite, and it is 
zero only for 6 = 3, 5, 7, . . . (or 6 = —2, —4, —6, . . . , but this case is excluded from our consideration). This 
observation, together with the behavior of (1 — x^Y''~^^^^ allows us to conclude that 



lim Q{x) = 



(X) if < 6 < 1, 
if 6 > 1, 



and the boundary condition (3.11 cannot be satisfied for any choice of C. We can therefore conclude that 
no solutions with 6 7^ 1, 2 exist. 



4. Analysis of the inviscid case: v — 



In this section we discuss the existence of solutions in case of zero viscosity. We show that the purely 
rotational trivial solution F = G = and = C^^ exists for all 6 > 0, but that no nontrivial solutions of 
the form (2.6 1 exist if 6 > 2. For 6 = 1, we present and discuss analytic solutions. For the cases < 6 < 1 



and 1 < 6 < 2, we present the governing equations that apparently admit nontrivial solutions. While we 
have not completed the existence and uniqueness analysis, we present our insights and numerical results in 
the case < 6 < 1, and we show that all potential solutions with 1 < 6 < 2 could not satisfy Rayleigh's 
circulation criterion and would thus be unstable with respect to axisymmetric perturbations. We again start 
by discussing the boundary conditions and then analyze the various cases that arise for various values of 6. 



4.1. Boundary conditions. Since slip is of no concern in the case of zero viscosity, we only need to address 
the no-source/sink requirements at the ground (a — tt/2) and at the center of the vortex {a = 0). From the 
analysis in the previous sections, it is clear that there will be no a priori restrictions on G or fl. Regarding 
restrictions on F, we have 



F ^ f = if 6 = 2, 
i^(0) = /(0)=0 and limF(a;)(l-x2)(i-'')/2 = lim/(a;) 







if 6 7^ 2. 



We still assume that near the vortex axis the azimuthal velocity, V0, behaves like C/r'', i.e., we still assume 
that (3.1) holds. We restate it here for completeness. 



lim n{x) = lim uj{x){l - x^)-'-^-''^/^ = G^ ^ 0. 



(4.1) 



4.2. Governing equations. Since the viscosity coefficient, is zero, (2.181 and (2.191 reduce to 

Csia) = 0, (4.2) 
C'i(a) + 26C2(a) =0, (4.3) 



which can be rewritten using (6.5) and (6.7) in the appendix. The third equation is (2.8) (restated later as 
(pl)). 



Remark 4.1. Note that the governing equations (2.3) and (2.4) in the case of zero viscosity are invariant 



under the transformation v i— )• — v, and so any obtained solution can be also "reversed" by changing its sign. 
In addition, we will see in some cases below that some of the equations are invariant under sign changes of 
some of the functions F, G, Q., etc., individually. 
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4.3. Case b — 2 (no nontrivial solutions). In this case, F = and (3.4| provides no information. Using 



(6.5) and (6.7 1, equations (4.2 1 and (4.3 1 reduce to 

G{x)Q{x) = 0, 

4.T 



(G^(x))' 



1 



(4.4) 
(4.5) 



Because of the boundary condition (4.1 1, we have that il 7^ in some interval (xq, 1) by continuity. Equation 
(|44| imphes that G = in {xq, 1), and ( [i^Sf reduces to (fi^)' = in (xq, 1). Thus = Ci^ in {xq, 1). This 
then implies that fl = Co, and G = in (0, 1). Thus we obtain the "trivial" solution 

F = G = 0, n = C^. (4.6) 



4.4. Case b ^ 2 (existence of the trivial solution). In this case, we can use relationship (3.4) between 



G{x) and F(x) and substitute it into (4.2) and (4.3). The resulting equations are given in the appendix in 
(6.13) and (6.15). Notice that we still have the trivial solution (4.6), since if = C^^, then (6.13) implies 
F{x) = cy/l — x^, and the initial condition F{0) = gives c = and F = 0; equation (3.4) then gives G = 0. 

There remains to be seen if there exist other, nontrivial solutions. We first address the simple case 6=1 
and then turn to the more complicated case b ^ 1,2. 



4.5. Case b - 



1 (existence of nontrivial solutions). If 6 = 

F{x)Q'{x) = 0, 
1 

2' 
2{FF"' 



1, then (6.13) and (6.15) reduce to 



{n'{x)y + ^(i-x'){F'{x)y" = o, 



(4.7) 
(4.8) 



where in we used the fact that (F^)'" = 2{FF"' + 3F'F"). 

We can consider two cases. Either F = 0, in which case (4.8) implies fl 



and from (3.4) we have G = 0. This case corresponds to the trivial solution m.G 



C^, (4.7) is trivially satisfied. 



In the second case, if F{xo) 7^ for some xq £ (0,1), then consider the largest interval (a, 6) C (0,1) 
containing xq such that F{x) 7^ in (a, 6) and F{a) = F{b) = 0. In (a, 6), ft has to be constant in view of 



(4.7) and F (x) = c{x — a)(b— x) with c > in view of (4.8). However, in this case all (one-sided) derivatives 



of F become infinite at a and b, and therefore the only possibility is that a = and 6 = 1, in which case we 

(1 



have n = Clj, 



Fix) = Ci^/x{l-x) with Gi ^ 0, and, from K^, G(x) = d 



In summary, when 6 = 1, we have a set of solutions of the form 

(1 - 2x)y/l + X 



2x)VT 

2y/^ 



F{x) = Ci y/x{l-x), G{x) = Gi - 



2^A 



for Gi e 



(4.9) 



which also includes the trivial solution (4.6) when Gi = 



We see that the solutions with Gi 7^ will have infinite fiow speeds both near the ground and near the 
vortex axis. The velocity becomes infinite near the ground in the radial direction (infiow for updraft solutions 
with Gi < and outfiow for downdraft solutions with Gi > 0) due to G having an asymptote a,t x = 0. Near 
the vortex axis {x — I) both fl and G have a finite limit and therefore the fiow speed becomes infinite due 
to the r term in the denominators of the velocity components (2.6). Both of these phenomena are observed 
in Fig. [2] in the middle plot. 



To assess how reasonable solution (4.9) is, we have implemented the iterative procedure described by 



Serrin |23| to compute solutions to (3.2) with small nonzero viscosity. In |23|, solutions depend on two 
parameters, k and P: viscosity is related to k via v — 2/k, so small values of viscosity correspond to large 
values of k. In Fig. [l] we present two solutions, one for k — 100 and P = (left) and one for k — 1000 and 
P = (right), and we compare them to (4.9) with Gi = G^j = 1. We observe very good agreement of the 
two solutions in the interval (0, 1) except near a; = 0, where the nonzero-viscosity solution exhibits a thin 
boundary layer due to the no-slip boundary conditions, while F' and G from (4.9) both tend to infinity as 



a; — >■ 0. Increasing k (i.e., decreasing the viscosity, v) results in shrinking the thickness of the boundary layer. 
This provides numerical evidence that downdraft solutions (Gi > 0) in (4.9) are limits of downdraft solutions 
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Figure 1. Comparison of the inviscid flow solution (4.9 1 with Ci = Ct^ — 1 with the 
solutions with small nonzero viscosity corresponding to fc = 100 and P = (left) and 
k — 1000 and P — (right) in (23]. The plot of F is shown in red, plot of G in green, and 



plot of fl in blue, with the plots for solution (4.9) dashed 






and isocurves for pressure (right) for the solution in (4.9 1 with Ci 
horizontal axes correspond to r = i? sin a and the vertical axes to 



Figure 2. Plots of the streamlines (left), the corresponding isocurves for speed (middle), 

4^2 and C„ = 1. The 
z = Rcosa. The values 
for speed range from 4 to 50 with a step of 1, increasing towards the vortex axis. The values 
for pressure range from —36.8 to 75.0 with a step of 2.3 with lowest values near the vortex 
axis. The straight red line corresponds to the level set p{R, a) = 0, which is the line 



cos a 



cl 



of (3.2 1 as — > 0. On the other hand, it follows from Serrin's results that updraft solutions (Ci < 0) in 
(4.9) cannot be limits as viscosity tends to zero of any of the solutions presented in [23]. This leaves open 
the question whether solutions other than those described by Serrin exist that tend to solution (14. 91) with 



Cl < as viscosity tends to zero. 

To visualize solution (4.9) in other ways, in Fig.[2|we show a streamlines plot, a contour plot of the speed, 
||v|| = vj^ + ''^a + ''-'e ! ^ contour plot of the pressure obtained from (2.16) with T — 0. The shown 
ranges are 0<r<l,0<2;<l with r = Rsina and z = Rcosa. We choose Ci = 4a/2 and C^^ = 1 
since the corresponding solution is also displayed later in red in Fig. |5] (as a limit of numerically computed 
solutions corresponding to b 1). 

We remark that the solution in ( |4.9[ ) is self-similar, and the self-similarity is clearly observed when one 
zooms out of the plots in Fig. [2] The zoomed out figures look identical to Fig. [2] with the contours only 
corresponding to different values for each level set. This is clear from the definition of the velocities in (2.6 1, 
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since the functions F, G, and ft only depend on a; = cos a. As an illustration of Cai's power law method, 
we computed the exponent for the velocity-radius power law from the data in the middle plot in Fig. [2] by 
computing the slope of the logarithm of the speed against the logarithm of the scale for several different 
pairs of points and obtained —1 for the slope. 

The pressure plot in Fig. [2]shows low values near the vortex axis and finite values along the ground, which 
increase as r — 0. In fact, from the solution (4.9l and the expression for pressure (2.16), one quickly obtains 
(taking T = 0) that 



piR,a) = - 



l-Cf{l 



l-Cf(l-cosa) 



SO we see that as one approaches the vortex axis, i.e., as a; — > 1, the pressure behaves like — It is also 
immediate to observe that the pressure is zero along the line 1 — Ci{l — x) — 0; this line is indicated by 
the bolder red line in Fig. [2] Notice that near the corner where the vortex axis meets the ground, all of the 
other level curves are tangent to this line, and there is a large pressure gradient as the pressure undergoes a 
sudden change from positive to negative values when crossing the red line and approaching the vortex axis. 
Clearly, such a behavior will only be observed if > 1, since only then does the level line p{R, a) = have 
positive slope. In Fig. [s] we show the cases that correspond to the line p{R,a) = having angles 7r/4, 0, 
and — 7r/4 with the horizontal, respectively. The corresponding values of are then 2 + \/2, 1, and 2 — y/2, 
respectively. 

The pressure plots in Fig. [2] and Fig. [3] also provide an interesting characterization of possible shapes of 
a visible tornado funnel. Since the funnel outline should approximately follow the isocurves, we see three 
distinct possible shapes: one that is conical near the ground (Fig. |2]and Fig. [s] (top row)); one that can be 
viewed as a degenerate, fully open cone, yet still with a single point touching the ground (Fig. [S] (middle 
row)); and one with the funnel having a nonzero width at the ground (Fig. [s] (bottom row)). Note that 
changing Ci while keeping Ci^ fixed changes the relative magnitudes of the azimuthal component of the 
velocity with respect to the non-azimuthal ones. We could view a large value of C^jCx as corresponding to 
large amount of swirl, and a small value corresponding to a small amount of swirl jsj T7j27 . In this sense, 
we can say that in Fig. |3] swirl increases from top to bottom, and wider funnels correspond to more swirl. 
We illustrate this behavior in Fig. |4] in which a few streamlines are shown for three cases, Ci = — 10 (left), 
Ci = -1, and Ci = -1/10 in (|4^. 



the lower-case functions and expressing g using (2.8 1, reduce to 

h 



4.6. Case 7^ 1,2 (general observations). In this case, the governing equations (4.2) and (4.3), using 

I, 
1 



2-6 



(4.10) 



(1-x^) 



2-6 



/7" + //" 



+ 4(1 - 6)//' + 2(1 - 6)(2 - fo) (/' + + 2(2 - = 0. (4.11) 



Remark 4.2. The last two terms containing uj in (4.11) can be written in terms of Q in the form 

= 2{2-b){i- x^)^-''nn', 



2(2 - b) 



ujuj' + (1 - 6) 



and hence we see from (4.11 1 that if F = f — 0, then we recover the trivial solution (14. 6|, F = G — and 



fl = Ccj. Similarly, if il = C^, then it follows from ( 4.10[ ) and the boundary condition F{0) = /(O) = that 
F = f = 0. We thus seek solutions with / ^O (or F ^ 0) and 17 ^ Ci^. 



It will be more convenient to work in terms of / and Q, so we rewrite equations (4.10) and (4.11) as 

b ' 



fn' 



1 



r + (2 - 6) 



X 



il-x') 



1-X2 
X 



f 



(4.12) 



-I- 4(1 - b)ff + 2(1 - 6)(2 - b)^^f'^ + 2(2 - b)(l - x'^)^^'' nn' = 0. (4.13) 

1 — x"^ 
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1 Illl^ / / 


















Figure 3. Plots of the streamlines (left), the corresponding isocurves for speed (middle), 



and isocurves for pressure (right) for the solution in (4.9) with C^j — 1 and Ci = \/2 + -\/2 
(top row), Ci = 1 (middle row), and Ci = ^2 — \/2 (bottom row). The values for speed 
range from 1.68 to 21.00 with a step of 0.42 (top), from 1.24 to 15.50 with a step of 0.31 
(middle), and from 1.30 to 13.78 with a step of 0.26 (bottom), increasing towards the vortex 
axis. The values for pressure range from —24.94 to 17.20 with a step of 0.86 (top), from 
—32.50 to (red line) with a step of 0.65 (middle), and from —0.7 to —35.0 with a step of 
0.7 (bottom), with lowest values near the vortex axis. 
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Figure 4. Illustration of the effect of the ratio Ct^/Ci on the swirl of solutions ( |4.9[ ). In all 
cases = 1, but Ci = —10 in the left plot, Ci = — 1 in the middle plot, and Ci — —1/10 
in the right plot. Four streamlines are plotted in each case, with initial points (20, 20, 0.02), 
(20,20,0.5), (20,20,2), and (20,20,5). 



Remark 4.3. Notice that both (4.12) and (4.131 are invariant under sign changes f i— — / and J7 i— — fl; 



as one consequence, we can assume that C^j > 0. 

On any interval where / ^ 0, we can solve ( |4.10[ ) for w and obtain 

uj{x) ^ c|/(x)p-'')/(2-'') for c e M. (4.14) 
Vice versa, on any interval where uj ^ (and in particular on some interval {xi,l) due to the boundary 



condition (4.1 1) we can solve for / and obtain 

fix) - c|c^(a;)p-f)/(i-'') = c(l - a;2)(2-f)/2|r!(x)|(2-'')/(i-'') 



for c G 



We thus have that 



O 



((1 - a;2)(2-b)/2) 



and 



F(x) = O (Vl 



1. 



(4.15) 



(4.16) 



4.7. Case 2 < b < oo (no nontrivial solutions). We now show that no nontrivial solutions of ( 4.12[ ) and 
(4.131 in the classical sense can exist for b > 2. First, from the definition of uj{x) = n{x){l — x^J^~^)/^ 
and the boundary condition (4.1 1 we see that uj{x) — > +oo as x ^ 1. Since in this case > 0, then, in 
view of (4.15); the only way to not violate the boundary condition f{x) as a; ^ 1 is if / is identically 



zero on some interval [xi,l). Since we are looking for solutions with / ^ 0, we can assume, without loss 
of generality, that / > on some interval {xq,xi). We now compare the limiting behavior of uj on either 
side of xi. On (xi, 1), where / = 0, we have from (4.13) that = C^^, so lo has a nonzero limit as a; — > 
from the right. On the other hand, on (a;o,a;i), where / > 0, we have (4.14), and hence lima;-;.^;^- a;(a;) = 0. 
Therefore lo cannot be continuous and there are no nontrivial solutions in the case 2 < 6 < cxi. 

4.8. Case 1 < 6 < 2 (behavior of potential solutions). In this case we have xE^ < 0. bmce again 



u){x) — +0O as a; — > 1, we can use (4.15 1 to express / in terms of lo and observe that this time the boundary 
condition f{x) — !• as a; — 1 is satisfied independently of the value of c in (4.15). Note that in this case lo 
has to be positive in (0, 1), since if a;(a;o) = for some xq e (0, 1), then, in view of (4.15), / would have an 



asymptote at xq and thus be discontinuous there. It follows from (4.15) that / cannot change sign in (0, 1) 
either, and (4.14) implies that ijj{x) — +cxi as a; — 0. Consequently, we also have ^{x) — > +oo as .a: — > 
and the azimuthal velocity becomes infinite near the ground. 
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We have been unable to find analytic expressions for such solutions and also encountered difficulties 
when approximating them numerically. However, in the next section, we show that such solutions would be 
unstable with respect to axisymmetric perturbations. 

4.9. Instability of potential solutions for 1 < 6 < 2. In this section we will assume < 6 < 2 and 
address the centrifugal stability of solutions to ( |4.12 ) and (4.13 1 with respect to axisymmetric perturbations. 
We will use Rayleigh's circulation criterion |7], which can be stated as the requirement that the Rayleigh 
discriminant $ is nonnegative, where 

1 d 



$(r) 



dr 



{rve) 



Substituting in the expression vg — Vl{x)/r^ and using the relationship between x and the cylindrical coor- 
dinates r and z, x = cos a — z/^/r"^ + z^, we obtain 



$(r) 



^2(1+6) 



Vl{x) [(1 - b)Vl{x) - x{l - x^)n'{x)\ 



(4.17) 



The particular case with b —1 gives $ = since from (4.9 1 we have VL = C^j- Also, the trivial solution with 
F = G = and 17 = C„ is clearly stable only for < 5 < 1. We will show that & = 1 is the largest value of b 
that permits stable nontrivial solutions. 

For 1 < 6 < 2, the stability requirement $(r) > can be replaced by an equivalent statement /^(a;)<i>(r) > 
since / 7^ in (0, 1). We can then use ( 4.12[ ) to get 

2{l-x^)fl^{x) 1-6 



f'ixMr) 



■fix) [{2~b)f{x)-xr{x)]. 



j.2{i+b) 2-b-' 

which is invariant under the sign change of /, so we can assume, without loss of generality, that / > in 
(0, 1). The stability requirement f^{x)^{r) > now implies (2 — b)f{x) — xf'{x) < 0, and, in particular, / 
is nondecreasing in (0, 1). However, this contradicts the assumptions / > in (0, 1) and lima;_j.i f{x) = 0, 
so no solutions corresponding to 1 < 6 < 2 can be stable with respect to axisymmetric perturbations. 



Lack of analytic solutions for 6=1 and constant nonzero 
We have not been able to find analytic 
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4.10. Case < 6 < 1 (numerical solutions) 

viscosity led to numerical approaches presented, e.g 
expressions for any nontrivial solutions in the case < 6 < 1 either, but we used a numerical approach to 
generate their approximations for various values of 6 between and 1. Once our solutions are computed, it 
can then be numerically or graphically verified that they satisfy the stability criterion $ > with $ given in 
All of our numerical solutions for < 6 < 1 were graphically checked to satisfy the correct inequality 



(4.171 



and thus were stable with respect to axisymmetric perturbations. 



We now describe our numerical approach to obtain approximations to solutions to (4.12) and (4.131. We 



first note that we can rescale the functions in consideration using the boundary condition (4.1) 



fix) = Cjix) 



and 



n{x) = c^n{x), 



so that the boundary condition (4.1) becomes 



lim D,{x) = 1. 



(4.18) 



Note that we can simply replace / and fl in (4.12) and (4.13) by / and fi, since the scaling constants cancel 
out. We will thus work with (4.12) and (4.13) in their original form and only replace (4.1) with (4.18), 



keeping in mind that any solutions will correspond to the rescaled functions. 

Disregarding the boundary condition on fl, it is clear that if the pair (/, il) solves (4.12) and (4.13), then 
so does any pair (±c/, ±££7) with c G M. In our numerical approach we will seek solutions with / > and 

> in (0, 1). Note that if / and fl satisfy (4.15), then equation ( 4.12[ ) will be automatically satisfied. It 
is not difficult to check that the pair 



/o(x) = 2(i-'')/2 {x{l~x)) 



(2-6)/2 



and 



flo{x) 



2x 



(l-6)/2 



(4.19) 
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satisfies (4.121 and the left-hand side of (4.131 evaluates to a well-behaved expression 

2 



-,1-b 



(2-6)(l-6) 



1 



that vanishes at both endpoints for all < 6 < 1 and converges uniformly to as 6 —> 1. The expressions 



in (4.19) can therefore serve as a basis for initial guesses in a numerical scheme. However, we observe that 



their derivatives behave singularly near the endpoints. To bypass this difficulty, we recall (4.16) and define 
a function 7(2;) by 

/(.T) = 7(x)(l-a;2)(2-'')/2, (4.20) 



so that, using (4.14) 



(4.21) 

We can then substitute (4.20) and (4.21) into (4.13). Since 7 is expected to have an infinite slope at 2; = 0, 



we also reformulate the newly obtained version of (4.13) in terms of the square of 7, 

p{x) = 7^(2;), 

and, after factoring out and discarding some positive terms, get the equation 



+ 2c2(l-6)p(3-2b)/(2-by 



(1 - x^){{l - x^)p"' - 2(4 - b)xp") - 2(2 + 6 - 3(2 - b)x')p' 

(1 - x'Wr - 2p ((1 - x2)p" - (2 - b)xp') 



(4.22) 



= 0. 



The boundary conditions arising immediately from those for / are 7(0) = and, in view of the asymptotic 
behavior of / given in (4.16 1, 7 having a finite limit as x ^ 1. Since the solutions can be rescaled as discussed 
above, we can assume 7(x) — >■ 1 as a; — > 1. We approximate the solution to (4.22) by discretizing it using 



a uniform mesh with stepsize h and solving the discretized system by Newton's method, using (4.19) and 
(4.20) to assemble an initial guess for p{x). The results presented in this section correspond to h — 10^^. 

While it is not clear to us whether solutions to ( |4.22 1 exist for any combination of the constants b and 
c, we have found that for a given value of b, increasing the value of c eventually creates instability in the 
numerical code, suggesting a potential restriction on (a combination of) these values. Note from, e.g., (4.20) 



and (4.21) that the constant c can be viewed as a scaling constant between and / (or, more generally. 



between the azimuthal component of the flow and the non-azimuthal ones); a large value of c corresponds 
to a relatively large azimuthal component of the velocity with respect to the other two components, while a 
small value of c corresponds to the azimuthal component being relatively small. In other words, increasing 
c can be viewed as increasing swirl in the flow. Notice the similarity to changing Ci in the case b — I above. 
In Fig. |5] we present computed solutions for c = 0.25 and b = 0.1, . . . , 0.9 with increments of 0.1. Notice 



that the results demonstrate continuous dependence on the parameter b, and the solution (4.9) for 6=1 can 
be viewed as their limit as 6 ^ 1. This solution, with constants C^^ — 1 and Ci = 4\/2, is plotted in Fig. [s] 
in red for comparison. We can observe that in most of the interval (0, 1) the magnitudes of F, 51, and G 
decrease as b decreases, although the results do not suggest that these functions would vanish if b approached 
0. Since the expression 1/r'' also decreases with decreasing 5 for < r < 1, we see that if c is fixed, then the 
flow speed decreases with decreasing b near the vortex axis. This is consistent with observations of Cai [2] 
and Wurman [28,|29) that larger values of b correspond to stronger storms. 

The azimuthal velocity exhibits an interesting feature for < 6 < 1. Notice in Fig. |5] that 57(0) = 
for every < 6 < 1, and therefore the azimuthal velocity vanishes at the ground. This behavior of 51 is 
not enforced by an a priori boundary condition, rather is it a consequence of the Euler equations and the 
boundary condition -F(O) — 0. It means that nontrivial solutions with < 6 < 1 exhibit purely radial inflow 
or outflow at the ground. 

To see how the choice of c affects the results, we also present results with a flxed value of b and varying 
values of c for which we were able to generate results. In Fig. |6] we present results with b = 0.6 and 
c = 0.1,. ..,1.0 with increments of 0.1. (The value b — 0.6 is chosen since it corresponds to the midpoint 
of the interval (—0.7,-0.5) found in (28].) For comparison, we also display the graphs corresponding to 
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Figure 5. Graphs of F{x) (left), ^l{x) (middle), and G{x) (right) as numerical solutions 
obtained from (|4.22[) with c = 0.25 for b = 0.1, 0.2, . . . , 0.9 (in blue). The functions F and 



f2 are increasing with increasing b at all x G (0, 1), while the magnitude of G increases in 



most of the interval (0, 1) as 6 increases. The red plot corresponds to the solution (4.9 1 with 
~ 1 and Ci = 4a/2, which can be viewed as a limiting case as 6 — 1. 




Figure 6. Graphs of F{x) (left), U{x) (middle), and G{x) (right) as numerical solutions 
obtained from ( |4.22[ ) with b = 0.6 for c = 0.1, 0.2, . . . , 1.0 (in blue). The functions F and 
n, are decreasing with increasing c at all x e (0, 1), while the magnitude of G decreases in 
most of the interval (0, 1) as c increases. The red plots correspond to the solutions with 
c = 0.25, which are also shown in Fig. [5] 



c = 0.25 in red; these same graphs are also shown in Fig.[5j We see that as c increases, the magnitudes of F, 
il, and G decrease in most of the interval (0, 1), although the change in fl is fairly small. This behavior is in 
agreement with the meaning of the constant c discussed earlier, i.e., that c reflects the relative importance 
of the azimuthal component of the velocity with respect to the other two components. 

To compare the flows, speeds, and pressure flelds, we present in Fig. [7] the analogs of Fig. [2] for the cases 
c = 0.25 and b — 0.8 and 0.2. The plots corresponding to the intermediate values of 6 showed continuous 
dependence on the parameter b and consequently we do not display them. We observe that while the 
streamlines remain relatively the same for various values of b, the speeds of the flow and the pressure fields 
exhibit discernible changes. 

Our numerical results in this section were obtained under the assumption that f2 does not change sign in 



(0, 1). We do not yet know whether our equations (4.121 and (4.131 admit solutions that change sign in the 
interval (0, 1). Such cases would be interesting to study, since, unlike in [23], equation (4.12) suggests that 



if / changes sign, so would f2 and vice versa. This could give rise to very interesting types of flows. 

5. Discussion, Conclusions, and Implications 

The phenomenon of vortex breakdown occurs where two conjugate flows, one super-critical and the other 
sub-critical, are matched. Vortex breakdown is a dynamic event, as opposed to our steady state model. A 
breakdown bubble is often found between the two flows. The super-critical flow, having a laminar appearance, 
erupts from the surface and does not allow waves to pass down from above to flU in the pressure deficit at 
its core. Surface friction inhibits the super-critical flow, and it is not in cyclostrophic balance, while the 
sub-critical flow above is in cyclostrophic balance and allows flow from above to fill the pressure deficit at 
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Figure 7. Plots of the streamlines (left), the corresponding isocurves for speed (middle), 
and isocurves for pressure (right) for the numerically computed solution with c ~ 0.25, 
b = 0.8 (top row) and c = 0.25, b = 0.2 (bottom row). The values for speed range from 
3.20 to 34.56 with a step of 0.64 (top), and from 1.7 to 6.6 with a step of 0.1 (bottom), 
both increasing towards the vortex axis. The values for pressure range from —24.0 to 54.4 
with a step of 1.6 (top), and from 1.82 to 14.56 with a step of 0.26 (bottom), with lowest 
values near the vortex axis. (The contours near the vortex axis for b — 0.2 are significantly 
affected by numerical errors in the contour plot routine.) The straight red lines correspond 
to the level sets p{R, a) — 0. 



its core. The region of the flow where the flow turns from horizontal into the vertical is called the corner 
flow region. Our model shares certain qualitative features with the corner flow region of the super-critical 
vortex below the breakdown bubble. We assume that the flow below the breakdown is quasi-steady and 
therefore use our model as an approximation. The vortices that we model have a singularity in the velocity 
field along the axis of rotation and along the surface near the base of the vortex axis. In atmospheric flows, 
the velocity field approaches the surface tangentially and is constrained to go to zero by the surface friction. 
In this work we provided numerical evidence that some of our solutions are viscosity solutions and hence 
the singular behavior near the ground can likely be removed by introducing a viscous boundary layer in 
which velocity tends to 0. Above this layer an inflow jet toward the axis can develop. Photogrammetry 
was used in [11' to document a horizontal inflow jet into a tornado with wind speeds of 55 m/s at treetop 
level in the dust band at distances up to 60 m from the dust column of the tornado. Along the vertical 
axis, numerical simulations suggest that a near-singular velocity fleld can develop [8][9|[l6]. There is also 
evidence that suggests that extremely intense vortices can develop into tornados by looking at subpaths 
(suction spots) within tornado paths in crops In |15j, the authors note that their solution was similar to 
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a similarity solution of Long 18 19 . Our solutions are motivated by Serrin's swirling vortex solutions 23 



which are conical similarity solutions, and Serrin's computations are analogous to Long's. Serrin's solutions, 
like ours, contain a singularity along the axis. Long's solution is not singular, however, he does not enforce 
a no-slip condition at the surface as Serrin does. Moreover, Shapiro and Markowski j24] comment that "the 
near-surface flow of Serrin's vortex beyond the core region may be a useful analog for the frictional boundary 
layer in the region of tornadoes beyond the radius of maximum wind." They suggest that a value of Serrin's 
result is the power drop-off in his solutions. We believe that current radar research in [2|l28 29 , giving the 



power-law drop for velocity where & ^ 1 (2 for vorticity) in tornados but varies, justifies our approach to 
modify Serrin's model. 

In this work, we focused on finding solutions of the form ( 2.6 1 to the Navier-Stokes and Euler equations in 
the upper half-space. We were motivated by the 1972 work of J. Serrin [23], in which he studies the viscous 
case with b = 1, and by later studies (2|[28 29 , in which it is suggested that the velocity may decay with 



different powers of the radial distance from the vortex axis than 1. We found that no solutions of the form 
(2.6) exist in the case of constant nonzero viscosity when 6 7^ 1, primarily due to the boundary conditions 
at the ground, i.e., the plane bounding our half-space. 

The situation is different in the case with zero viscosity and Euler equations. In this case the boundary 
conditions are relaxed since slip is allowed. Assuming that the azimuthal velocity behaves like C^/r'' with 
near the vortex axis, we were able to show that the trivial solution = and F = G = works 
for any & > 0, although it is stable with respect to axisymmetric perturbations only if < 6 < 1. 

Nontrivial solutions are harder to find. We showed that if 6 > 2, no nontrivial solutions exist. We 
also showed that any potential solutions for 1 < < 2 would be unstable with respect to axisymmetric 
perturbations. The case with 6=1 was fully analytically resolved and its solution is given in (4.9 1. We 
discussed its characteristics and showed that the downdraft (Ci > 0) solution in (4.9) can be viewed as a 
limit as viscosity goes to zero of the downdraft solutions found in [23] . 

The case with < & < 1 proved analytically difficult, and we only presented some numerical results 
indicating the existence of solutions that have similar characteristics to those with 6 = 1. This case is most 
interesting, as it allows the coefficient 6 to fall into the ranges discussed by Cai and Wurman et al. |2, 28,29.. 
We were able to numerically find solutions with F and ft that do not change sign. It would be interesting 
to see whether solutions with sign changes are possible. 

Finally, we remark that in the case of a turbulent motion the functions of the form (2.6 1 correspond to the 
mean velocity, and thus v could play a role of eddy viscosity to maintain the delicate balance between the 
mean and turbulent components of the flow. In recent studies, testing eddy viscosity assumptions with direct 
numerical simulations showed varied success [S], but has remained an important tool for understanding the 
connection between the scale of the model and dissipation of energy 13 . Following [22], Serrin suggests an 
experimentally motivated relationship v (T||v||r, where cr is a small dimensionless constant 23 . Taking into 



account (2.7), we thus obtain u ■ 



r, where cr is a small dimensionless constant 23 
with the simplest case being r 



const. Such an assumption on 



viscosity then leads to a modification of equations (2.18) and (2.19) that will be investigated in the future. 



6. Appendix 



6.1. Navier Stokes equations in spherical coordinates and incompressibility. The three compo- 



components (2.5) have the form 



nents of the Navier-Stokes equations (2.3) expressed in spherical coordinates and in terms of the velocity 



dvT. 



Va dVf 



+ 



Vg dvji 



dR R da i?sina dO 



R 
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dvg Vadva vg dvg VRVg - vl cot a 
dR R da Rsina 89 R 
I dp V 



Rda i?2 



R 

OR \ dR J sin a da 

1 d^Va ndVR 



da 



sin^ a dO"^ 



da 



sm a 



Va + 2 COS a 



dv0 
1)0 



dv0 Va dvg 



Vg dvg VRVg + VgVg COt a 



dR R da i?sinQ! dd 
1 dp v 



Rsina de R^ 



R 

dR \ dR 



1 d f . dvg 
sm a da \ oa 



1 d^ve , 1 
sin2 a 



sin^ a 



r, • dvjj dva 



Similarly, the continuity equation (2.4 1 has the form 

WdR^^ ''^^ + i?sin^ 



d . .dvg 
-Ksma) + — 



= 0. 



(6.2) 



(6.3) 



(6.4) 



6.2. Expressions Ci and Di for general 6 > 0; the continuity equation. Substituting variables (2.6 1 
into the Navier-Stokes equations (6.1|-(6.3) and comparing with the forms in (2.12)-(2.14), we obtain the 
expressions for Ci{a) and Di(a) (we omit the argument x = cos a in the functions F, G, fi, /, g, and w, 
and, for example, we write G' instead of G'(cosa)) 

F"^ + bG'^ + + hFGcoia^ FG' sin a 



Ci{a) - 
C2{a) = 
C3(a) = 



2b 



(sin a) 

-{bF^ + n^)cota+{l-b)FG- FF'sina 
(sinQ;)^^ ' 
(1 - b)Gn + F[{1- b)fl cot a-Q' sin a] 



(sin a) 



2b 



and 



Diia) 



D2{a) 



Dsia) 



G" sin^ a - 2(1 - fe)G' cos a - (1 - 62 - cos 2a)G csc^ a - 2(1 - 6)F cot a + 2i^' sin a 

(sina)^ 

F" sin^ a - 2(1 - 6)i^' cos a - (1 - b^)Fcac'^ a - 2& G cot a - 2G' sin a 

(sina)'' ' 
sin^ a - 2(1 - b)fl' cos a - (1 - b'^)n csc^ a 



In terms of x, and the functions F, G, il, and /, g, w, these expressions can be written as 



Gi(x)-(l-:r2)-b 
= (l-x2)-i 

C2ix) = {1 - x'r" 



-F^ -bG^ -n'^ - F ( Vl - x^ G' + b 



-,f-bc/-oj^-f[Vl~x'9' 



Vl-x 

X 



G 



x^ 



9 



{bF^ + f]^) - F Vl - a;2 F' - (1 - b)G 



-^^^if + ^')-f {Vl^f - (1 - %) 
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C:i{x) = [1 - x') 



2\-b 



(1 - h)GVL - F ( - a;2 fl' - {1 - b) 



and 

^{l-x')-'^' 



(1 - x^)-^ \{1 - b)guj - v/l-a;2 fuj' 



2(1 -x^)- 62 



(6.5) 



(1 - x^)G" - 2(1 - b)xG' - " ' ^ G - 2(1 - h) , F + 2^/l^ F' 

(1 - x^)g" + ~ (2 - 6)(1 + b)]g + 2^/l~x^f 



D^ix) = (l-x2)-^/2 

= (l-x2)-^/2 



(1 - x^)F" -2{\-b)xF' - - — F -2b 



1-62 
1 -a?2 



G-2^/\-x^^ G' 



(1 - x')f" - 6(1 - 6)/ - 2 



2\/\- x^ g 



2 „' 



1-62 ■ 

(1 - x'^)Vl" - 2(1 - b)x Q! - ^ n 

1 — 



2^3(2:) = (l-x2)-''/2 

= (1 - x2)-''/2 [(1 - V - 6(1 - 6) uj 
The expressions Ci + 26 C2 and £)i + (1 + 6)Z?2 are then 



(6.6) 



C1 + 26C2 = (l-x") 



2\-b 



26 



((1 - 6)f2 + 6 G^) + 2v/l - ((1 - + 6 GG' + VlVl') 

FG + b X F'G - {I - 3b)x FG' 



+ (l-.^)(^'G' + ^G") + 5^-^^-2(l-^^)^^ 



1-X2 



= (l-x^ 



2x 



((1 - h)f + bg^ + (1 - 6)c^2) + 2v/l-a;2 ((1 - 6)//' + 6.95' + ljuj') 



(6.7) 



2..,, , , , l + 2.x2 + 26(l-6)(l-a;2) 



+ (l-:r2)(/V + /g") + 



l-a;2 



fg + xfg + 2xfg' 



and 



Di + il + b)D2 = il-x 



2\-{2+b)/2 



(1 - 6)(1 - 6^ - 26(1 - x'^))F - 6^(4 + 6 - 2^^) 



G 



+ 2(1 - 6)2^(1 - x^)F' + (2 - 46 - 6^ - 2(1 - 36 + b^)x^)Vl - x^ G' 
-{l-xY 



(1-^2) 



2N-3/2 



(1 - b)F" - (4 - 36) ^ G" + G" 

V 1 — x2 



6(1 - 62)(1 - a;2)/ + (-3 - 6(1 + 6)(1 - x^)) 



(1 + 6(1 + 6)(1 -a;2))Vl-x2 g' 



2\2 



(1-^2) 



(1 - b)f" 



Finally, substituting ( |2.6| and ( |2.7| into the continuity equation ( |6.4| , we get 

X 



i?-(i+b)(l_a;2)-V2 



(2 - 6)G - (^v/r^F' - (1 - h)-jJ= F^ 



0, 



_ 2:2)^1/2 [(2 _ b)g _ - a;2 /'I = 



(6.8) 
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6.3. Expressions C,; and Di for 6 7^ 2. In this case, we can use equations (2.8) and eliminate G{x) and 
g{x) from the expressions in the previous section. We have 

'1-b 



= (1-X2)-V2 



2-6 
1 - b 



n{x) { F'{x) + ^— F{x) ] - F{x)n'{x) 



2-6 



fixMx) - fixWix) 



D3 = (1 - [(1 - x2)f7"(a;) - 2(1 - 6)x(l - x'^)n'{x) - (1 - 62)f](a:)] 

= (1 - _ a;2)cj"(.T) - 6(1 - 6)w(x)] , 



C1 + 26C2 



(^_^2)l/2-b 

(2-6) 



(1 - ( '^±^F'{x)F"{x) + F{x)F"'{x)^ + 2(2 - b)n{x)n\x) 
^#±g^-(x) + ^+f + f-V (.)F-(.) 
+ (2 + 6)a;(i^'(a;))2 + (4 - 6)a;i^(a;)F"(a;) 



1- 6 

2- 6 



(l-a;2)-i/2 



2-6 



(1 - (|^/'(:^)/"(^) + /(2:)/"'(x)) + 2(2 - bMxWix) 



and 



Z)i + (1 + b)D2 



+ 2(1-6) 

(1 - a;2)-2-6/2 



(2 - b)Y^ (/'(^) + + 2/(2:)/'(x) 



2-6 



(l-x2) 



2N-1/2 



2-6 



(1 - xYF^^H^) - 4(2 - 6).i;(l - xyF"'{x) 

- 2(1 - 6)(3 + 6 - 2(3 - 6)a;2)(l - x^fF"{x) 

- 46(1 - 6)(2 + 6 - x^)x{l - x^)F\x) 

- (1 - 6)(3 - 6(1 - 6 - 6^ - 4(3 + b)x^ + 4x^))F{x) 
{l~x'rf<^^\x)-4x(l-x')r{x) 

- 26(1 - 6)(1 - x^)r{x) + 6(1 - 62)(2 - b)f{x) 



(6.9) 



(6.10) 



(6.11) 



(6.12) 



6.4. Equations ( |3.3| for 6 7^ 2. When 6 7^ 2, we can use expressions ( |2.8| and substitute them into 



equations (3.3 1 to get, in terms of F and fl, 



Fn' = 



1 - 6 



2-6 



l-x2 



2(2 - 6)r2r2' + (1 - 

1 - 6 



(1 - a;2)2r2" - 2(1 - 6)a;(l - x^)n' - (1 - 6^)17 = 0, 
^ ^ ^ F'F" + FF'" 



= 2 



2-6 



2-6 

a;(l + 6a;2) 3 6 + (2 + 36)a:2 



(1 -a;2)2 



F' 



1-x^ 



FF' + (2 + 6).t(F')^ + (4 - b)xF{x)F" 



(1 - - 4(2 - 6)a;(l - x^fF'" - (1 - 6) [2(3 + 6 - 2(3 - b)x'^){l - x'^ fF" 

+ 46a;(2 + 6 - x^){l - x^)F' + (3 - 6 + 6^ + 6^ + 46(3 + h)x^ - 4bx^)F 



= 0, 



(6.13) 
(6.14) 

(6.15) 
(6.16) 
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or, in terms of / and u>, 

(1 - x^)uj"{x) - 6(1 - b)uj{x) = 0, 
+ 4(1 - b)ff' + 2(1 - 6) (2 - b)^-^ (f + Lu^) + 2(2 - = 0, 

(1 - a;2)V^*' - 4a;(l - x'^)f" - 26(1 - 6)(1 - a;^)/" + 6(1 - b^){2 - h)f = 0. 
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